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Abstract. 

In this paper we report some important results that help in analizing the border 
collision bifurcations that occur in n-dimensional discontinuous maps. For this purpose, 
we use the piecewise linear approximation in the neighborhood of the plane of 
discontinuity. Earlier, Feigin had made a similar analysis for general n-dimensional 
piecewise smooth continuous maps. Proceeding along similar lines, we obtain the 
general conditions of existence of period- 1 and period-2 fixed points before and after 
a border collision bifurcation. The application of the method is then illustrated using 
a specific example of a two-dimensional discontinuous map. 
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1. Introduction 

In recent times a significant amount of research effort has been directed towards 
nonsmooth bifurcations because a large number of physical and engineering systems 
have been found to exhibit this class of bifurcation phenomena pi EJ Ej . Such physical 
systems are generally given by two or more sets of differential equations, and the system 
switches from one set to another when some switching conditions are satisfied. It has 
been shown that such systems give rise to piecewise smooth maps, where the discrete- 
time state space is divided into two or more compartments with different functional 
forms of the map, separated by borderlines [HE]- A new type of bifurcation, called border 
collision bifurcation occurs when a fixed point collides with a borderline, resulting in an 
abrupt change in the Jacobian matrix. Most past studies on border collision bifurcations 
have concentrated on piecewise smooth maps that are continuous across the borderlines 

B El U EDI EEH- 

However, it has also been found that some switching dynamical systems, for 
example, the Colpitts oscillator [12] . dc-dc converters [T3], thyristor controlled reactor 
circuits [14] . sigma-delta modulators 15J and many other electronic circuits [16J give 
rise to discontinuous maps. Discontinuity in the map arises if, in the n dimensional 
Poincare section there exists an (n — 1) dimensional manifold such that infinitesimally 
close points at the two sides of the manifold map to points that are far apart. In order to 
explain the bifurcation phenomena in such systems, it is necessary to have a bifurcation 
theory for discontinuous maps. 

The bifurcation theory for ID and 2D continuous piecewise smooth maps has been 
developed in [HI EJ. In another important line of work, Feigin obtained some results 
on the existence of periodic orbits in n- dimensional piecewise smooth continuous maps, 
which were brought to the notice of the English speaking audience in The theory 
for discontinuous maps is in the preliminary stage of development, with some progress 
reported for one-dimensional maps in [T7] ITS] . 

In the present paper, we follow Feigin's approach in case of discontinuous maps, 
and present new results that are applicable to the general n- dimensional case. In 
particular, we present the conditions of existence of the period- 1 and period-2 fixed 
points which form the basis of the observable bifurcation phenomena in such maps. For 
this purpose, we use a piecewise linear normal form where the discontinuity is defined 
by two parameters. 

Finally, the application of these results is illustrated by deducing the bifurcation 
behavior of a two-dimensional discontinuous map, and validating the results by 
simulation. 
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2. The Normal Form of the Discontinuous Map 



We consider discontinuous maps of the form 

fX fc 



X 



(fc+1) 




for XeP- 
for X G P 



where the discrete-time phase space is divided into two regions P~ and P + , and II ~ and 
Il + represent the functional forms of the map in the two regions respectively. 

Since our points of interest are the occurrence of bifurcations with the traversal 
of a fixed point across the border, we can study the same by applying local linear 
approximation in neighborhood of the discontinuity: 



(fch 



where 



and 



x(fc+ i) = 



Ai 



AiXW + Cm, ifx[ k) < 



A 2 X« + C//2, Hx[ k) > 



(1) 



<9ir 



ax 



A 9 



x=o 



<9X 



x=o 



<9IT 



(ii=0 



d/i 2 



M2=0 



Here X^ fc+1 ^ is (k + l)th iteration of n x 1 matrix X = (xi, x 2 , 



Ai and A 2 are 



n x n matrices, C = (ci, c 2 , 



is n x 1 matrix and \L 2 are parameters. 



We can apply a suitable coordinate transformation to obtain the matrices in the 



normal form, where the matrices Ai = [aj£ ] and A 2 = [af)] are such that = 



if j ^ 1, i.e., except the first columns of the matrices all other columns are same 
|| . In the earlier work on piecewise smooth continuous maps |0j this condition was 
required to ensure continuity of the map. In the present work, though we are considering 
discontinuous maps, we retain this property of the matrices Ai and A 2 and move the 
discontinuity to the parameters Hi and /i 2 without loss of generality. In [T7j the length 
of the discontinuity was considered to be constant. In this paper, in order to allow a 
variable length of the discontinuity, we have used two parameters n\ and \i 2 . We have 
assumed C to be same on both sides, which physically means that the two parameters 
have proportional effect on the two sides of the border. 



2.1. Existence of period- 1 fixed points 

Let us consider two fixed points X^ = {x* lL , x^ L , ■ ■ ■ , x* nL ) T and = 
(x* R , x 2i? , . . . , x* nR ) T of the sub-mappings II - and Il + respectively. That is, 

X£ = AiX£ + C^i for x* 1L <0, (2) 

X* R = A 2 X* R + C^ 2 for x* 1R >0. (3) 

|| The derivation of the normal form in one-dimension can be found in 8 , and that in two-dimension 
can be found in ^3 
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where x\ L and x* lR are the first components of and respectively. 

We define p*(X) = |AI — AjJ and p**(X) = |AI — A 2 |, where I is the n x n identity 
matrix ^f. 

Assuming (I — A x ) and (I — A 2 ) to be invertible, we find 
XI = (I - AO- 1 ^ = s^ffiCfr, 

Now the kth element of these matrices can be written as 

Kl Kr 
x kL = ^(I)^ 1 ' x hR = ^(i)^ 2 ( 4 ) 

where b* kL = [(adj(I — Ai))C]fe is the kth element of the column matrix (adj(I — Ai))C 
and b* kR = [(adj(I — A 2 ))C]fc is the kth element of the column matrix (adj(I — A 2 ))C. 
Now, since in the matrices A x and A 2 , = a\ 2 if j 7^ 1, the cofactor of first column 
of (I — Ai) and (I — A 2 ) are same. Hence 



i=l i=l 

where and , are the elements of the first rows of the matrices (adj(I — Ax)) and 
(adj(I — A 2 )) respectively and q are the elements of the column matrix C. 
Using the above in (@J) we obtain 

b* h* 
XlL ~ ^(l)^ 1 ' XlR ~ p**^)^ 2 { > 

We require x* XL < for X^ to exist and x* 1R > for X.* R to exist. With the help 
of above derivations we can now attempt a classification of the possible behavior at the 
border collision bifurcation. 

Equation (JSJ) points to two possibilities: 
Case 1: p*(l)p**(l) > 0. This will be true if both p*(l) and p**(l) are positive or if 
both are negative. There are two sub-cases. 

Case la: If /zi/i 2 > then only one fixed point will exist. 

Case lb: If /ii/i 2 < then either two fixed points will exist or no fixed point will exist. 

Case 2: p*(l)p**(l) < 0. This will be true if p*(l) > and p**(l) < 0, or if p*(l) < 
and p**(l) > 0. Here also there are two sub-cases. 

Case 2a: If /^i/x 2 > then either two fixed points will exist or there will not be any 
fixed point. 

Case 2b: If \x\\i2 < then only one fixed point will exist. 

Table [T] and Table El summarize these results for b\ > and b* < respectively. 
From these tables the bifurcation sequences can be inferred for a given set of parameters, 
as /ii and /i 2 are varied. 

% Note that Feigin JI] defined these in the opposite way, as p*(A) = |Ai — AI| and p**(\) = IA2 — AI|. 
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Table 1. The conditions of existence of fixed points when b\ > 0. 



Conditions on the parameters 


Existence of fixed points 


p*(l)p**(l) 




p**(l) 


Mi > 

M2 < 


Mi > 
M2 > 


Mi < 
M2 > 


Mi < 
M2 < 


> 


< 


< 


X* L kX* R 


xi 


None 


xi 


< 


> 


< 


XL 


None 


xi 


xi & xi 


> 


> 


> 


None 


xi 


xi &xi 


xi 


< 


< 


> 


X* L 


xi & xi 


xi 


None 


Table 2. The conditions of existence of fixed points when b\ < 0. 


Conditions on the parameters 


Existence of fixed points 


p*(l)p**(l) 




p**(l) 


Mi > 
M2 < 


Mi > 
M2 > 


Mi < 
M2 > 


Mi < 
M2 < 


> 


< 


< 


None 


xi 


xi &xi 


xi 


< 


> 


< 


X* L 


xi&xi 


xi 


None 


> 


> 


> 


x* L &x* R 


xi 


None 


xi 


< 


< 


> 


xi 


None 


xi 


xi &xi 



3. Conditions for the existence of period-2 orbit 

Suppose that a new two period solution originates at a border collision point. Let the two 



**\T 



points of the period-2 orbit be M* = (m*, m*,, . . . , m* n ) T and M** = (ml*, m* 2 * , . . . , m**) 
lying on the half planes L and R respectively. 

Now we derive the condition for the existence of the period-2 orbit. Using the local 
map we can write 



f(M*) = M** and /(M**) 



M* 



i.e., 



M* = A 2 M** + C// 2 , 
M** = A X M* + Cmi- 
Now subtracting (JHJ) from (J7J) we get 

AM = AxM* - A 2 M** + C(mi - 

where AM = M**-M* = (m\* -m* 1: m* 2 * -m* 2 , . . . , m**-m*) T = (5m 1 ,5m 2 , 
Now as ai} = af) = a^- if j ^ 1, writing (jHJ) in scalar form we obtain: 



Kl 



5mt 



,(2)„ 



a kl m l ~ a kl m i 



+ c k (ni - M2) + 2J a fcj( m j - 

3=2 

for k — 1, 2, . . . , n. 

Equation can be rewritten by adding and subtracting a^ml* as 



(6) 
(7) 

(8) 
5m n ) T . 

(9) 



8mi 



'ki 



(2)\ 
a kl' W 



T + c fc(^i - ^2) + 5^ a kj( m ) 
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Similarly by adding and subtracting a^}m\ in equation (JHJ) we get 

n 

Sm k = (4i - 4? ) m i + c *(^i - ^2) + 2 a kj( m j - m */y 

The above two equations can be written in vector form as 

AM + A X AM = [Ax - A 2 ]i m\* + C(/xi - // 2 ) (10) 
AM + A 2 AM = [Ai — A 2 ]i m\ + C(^i - fi 2 ) (11) 

where [Ai — A 2 ]i indicate the first column of [Ai — A 2 ]. 

Assuming (I + A 2 ) to be invertible, from ()11|) we obtain: 

AM = ^)^ + ^T)^-^ (12) 

where B* = (-l) n (adj(I + A 2 ))[Ai - A 2 ]i and D* = (-l) n (adj(I + A 2 ))C. Similarly 
assuming (I + Ai) to be invertible, from (fTTH) we obtain 

AM = ^zyyK* + - *0 ( 13 ) 

where B** = (-l) n (adj(I + Ai))[Aj - A 2 ]i and D** = (-l) n (adj(I + Ai))C. 
Now writing (12) and (13) in scalar form 

6mk = p^T) ml + ^Fi) ~ ^ (14) 

6nik = ^TT m i* + -^TT(^i-^2). (15) 
p y — 1) p (— lj 

Since ajy = ajy = for j ^ 1 and & = 1,2, ... ,n, we have &** = b\ = b\ and 
rf** = d\ = d\. Hence from (14) and (15) 

h d x h m dx 



p**(-l) P*(-l 
where /x = (/ii — /i 2 ). From (14) 

/)* d* 

m k = m h + ^rry^i + ^piy ^ (17) 

substituting (17) into (7), we obtain (in scalar form) 

* (1) * * 

m k = }_^ a kj m j - p ^_^ m i - ^ZTJ ^ + c fc /n 

or equivalently using matrices: 

M* = A X M* + Cah - ^yy^o (18) 
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where 



( b\ ■ ■ ■ \ 



A x = Ax 



1 



p* 



b\ ••• 



b* ••• 



\ o 

\ n 



Assuming (I — Aj) to be invertible, we can now solve M*. From (18) 



M* 



Now 



where 



II- Ail 



1)^° 

|Ip**(-l) - Aip**(-1) + L 



d 



1) 



/ 6* 
6* 



1) 



0\ 




V o • ■ • o J 



and d is some constant. 

Hence from f!2()l) and (|21j) we obtain 

adj^-AO^-l) 



M* 



d 



C/i! 



D* 



p**(— 1) 



/to 



where 



~ adjQt-AQC ~ adjCI-AOD* 

Ki = ; and K2 = — 



(19) 



(20) 



(21) 



(22) 



d d 
Writing (J2~2")) in scalar form 

m* = kxnip**(-l) - k 2 {ni - /i 2 ), (23) 

and substituting this into (16) we obtain 

ml* = - fc 3 (/ii - fi 2 )- (24) 

Hence in two parameter n-dimensional discontinuous map (JIJ, a period-2 orbit will 
occur only when the following two conditions are simultaneously satisfied: 

fcip**(-l)/zi - hi/it - fr) < 0, (25) 
fciP*(-l)//i-fc 3 (/ii-// 2 ) >0. (26) 
where fci and &2 are respectively the first components of Ki and K2, and 



fa = h 



P *(-r 



fin 



: p**(-i) 1 &! y ' 

The derivation of the conditions for the existence of orbits of periodicity greater 
than 2 in the n-dimensional discontinuous map is very cumbersome. This is still an 
open problem. 
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Note that (J25|) and (|26|) give the conditions of existence of period-2 orbit. In 
addition, one has to consider the stability of the orbit in order to determine whether 
the period-2 orbit will actually be observable in the asymptotically stable behavior of 
the system. 



4. Example: A Two Dimensional Piecewise Linear Discontinuous Map in 
Normal Form 

The normal form for the two dimensional discontinuous map is 

A 1 



y (n+l) 



< 



A 9 




Cnx, if < 



(2f 



where 



Ai 




A 9 



+ C/x 2 , if x {n) > 



C 





Tl is the trace and 5l is the determinant of the Jacobian matrix Ai of the system at a 
fixed point in L := {(x, y) G 3ft 2 : x < 0} and tr is the trace and 8r is the determinant 
of the Jacobian matrix A 2 of the system at a fixed point in R := {(x, y) G 9ft 2 : x > 0}. 
The normal form (J28|) follows from that derived for piecewise smooth continuous maps 
in [3110], with discontinuity added through the parameters \i\ and 

In case of the map (|28J). the terms p*(l), p**(l), p*(— 1) and p**(— 1) are given by 



P*(l) 



P 



P 



II- Ail 



II- A, 



T — Ai| 



I- A, 



l-r L -1 
h 1 

1-TR -1 

5 R 1 

-1-T£ 

5l 



1-TR -1 
S R "I 



(l-r L + 5 L ), 

(1-Tfl + ^). 
= (1 + T L + 5 L ), 
= (1+T R + S R ). 



4-1. Inference about period-1 fixed points 



Now if we choose p*(l) > 0andp**(l) > 0, i.e., p*(l)p**(l) > 0, application of the results 
of n-dimensional discontinuous map (for n = 2) will yield the following inferences. 
If we keep fi 2 constant and vary fj,i, then 

• if > then only one period-1 fixed point will exist, 

• if /ii//2 < then either two period-1 fixed points will exist or no period-1 fixed 
point will exist. 
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We now show, using the specific 2D map (J28j) that the above results are true. 



4-2. Proof of the existence ofX.* L and fixed points 



If a fixed point X£ 
equations 



x* Ll y* L ) T exists, it is given by the solution of the system of 



X i 




Now the fixed point will exist if x* L < 0, i.e., 
^ = ^<0 

l-T L + S L P *(l) 

This is true if 

{tl>1 + 5l or equivalently p*(l)<0} for /ii>0, 
{tl<1 + 5l or equivalently p*(l)>0} for fii<0. 
The above conditions (j3~Tf and (J32|) are depicted graphically in Fig. QJa) 



(29) 

(30) 

(31) 
(32) 



Sl 



8 R 



1.5 
1 

0.5 





-0.5 
-1 

-1.5 



-1.5 -1 -0.5 



(a) 

m<o 




i/ 


/ Ui>0 



1.5 
1 

0.5 



0.5 1 1.5 



-1 
-1.5 



(b) 
|^2>0 






/ l^2<0 



-1.5-1 -0.5 



0.5 1 1.5 



Figure 1. (a) Regions of existence of the X£ fixed point as /ii is varied through 0. 
The fixed point exists in the shaded region if \i\ < 0, and in the unshaded region if 
Hi > 0. (b) Regions of existence of the X^. fixed point as /i2 is varied through 0. The 
fixed point exists in the shaded region if [i^ < 0, and in the unshaded region if fi2 > 0. 

In a similar manner it can be shown that the fixed point XJj will exist if 

; *TT X =^7TT>0 (33) 
l-T R + 5 R p**{l) 

This is true if 

{t r >1 + 5 r or equivalently p**(l)<0} for ^ 2 <0, (34) 
{t r <1 + 5 r or equivalently p**(l)>0} for ^ 2 >0. (35) 

The above conditions (}34"j) and (}35|) are depicted graphically in Fig. HJb) . 
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The conditions (|31|). ()34|). and (J35|) divide the (t l ,t r ) parameter space into 

four regions R\, R2, R3 and R4 as shown in Fig. |21 These regions correspond to the four 
rows in Table d 



* 1+5 L 



R 2 




Ri 








R 3 

1 


(0,0) 


R 4 



Figure 2. Partitioning of the parameter space where different bifurcations should be 
observed. 

In the region R ly since r L > 1 + 5l and r R > 1 + 5r, the period one fixed points 
and "K R are not stable. Therefore, even though they exist, the trajectory diverges 
to infinity. 

0.15 
0.1 

0.05 

x 



-0.05 
-0.1 

-0.1 -0.05 0.05 0.1 

^1 

Figure 3. A one-parameter bifurcation diagram with /ii as the bifurcation parameter 
and fi2 = —0.025. The other parameters tl — 0.3, 6l = 0.3, tr = 1.5 and 5r = 0.3 
are in R2. 

In i? 2 , suppose we keep ^2 fixed at a negative value, and vary Tabled shows 
that the period-1 fixed point will exist so long as \i\ < 0, and the fixed point 
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^1 



Figure 4. A one-parameter bifurcation diagram with /ii as the bifurcation parameter 
and 1x2 = 0.025 . The other parameters Tl = 0.3, Sl = 0.3, tr = 1.5 and Sr = 0.3 are 
in R 2 . 

X^j will exist for the whole range of fj,\. Now, if we choose the parameters such that 
— (1 + Sl) < tl < 1 + 5l and Tr > 1 + 5r, then is stable and X.* R is unstable. The 
representative bifurcation diagram for such a situation is shown in Fig. El which shows 
the period- 1 orbit for n\ < which disappears for /i 2 > 0. If we take /i2 > and vary 
fix, then we infer from Table^that X^ will exist (and is stable) for fi\ < 0, but no fixed 
point will exist for li± > 0. This is what we see in the bifurcation diagram of Fig. EJ 




-0.2 -0.1 0.1 0.2 0.3 0.4 0.5 



^1 

Figure 5. A one-parameter bifurcation diagram with /ii as the bifurcation parameter 
and fi2 = —0.025 . The other parameter T£ = 0.7, Sl = 0.3, tr = 0.3 and 5r = 0.3 
are in R3. 

In the region i? 3 , suppose we fix jj, 2 < and vary Hi. Table [T] shows that for ^ < 0, 
X^ will exist but X.* R will not. For fii > none of the fixed points will exist. Now, if 
the parameters are chosen such that —(1 + 5 L ) < t l < 1 + then a stable period-1 
behavior will be observed so long as ll\ < 0. This is shown in Fig. El On the other hand, 
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0.04 



0.02 



-0.02 



-0.04 




Figure 6. A one-parameter bifurcation diagram with /ii as the bifurcation parameter 
and fi2 — 0.025 . The other parameters T£ = 0.7, Sl = 0.3, tr = 0.3 and Sr = 0.3 are 
in R 3 . 

if fi 2 is chosen to be positive, for /x x < both the fixed points will exist, and for \i\ > 
X|j will exist but will not. If the parameters are such that —(1 + Sr) < tr < 1 + Sr, 
then XJj will occur for the whole parameter range, while X^ will occur only for lli < 0. 
This is shown in the simulated bifurcation diagram of Fig. |H1 



0.35 
0.3 

0.25 
0.2 

0.15 
0.] 

0.05 


-0.05 



-0.1 



R 4 ( ^2 < ) 




-0.6 



-0.4 



-0.2 



0.2 



0.4 



0.6 



Figure 7. A one-parameter bifurcation diagram with /ii as the bifurcation parameter 
and fi2 = —0.025 . The other parameters T£ = 1.5, Sl — 0.3, tr = 0.3 and Sr = 0.3 
are in R4. 

Similarly in R4, suppose we keep \x-i fixed at a negative value, and vary Tabled 
shows that the none of the fixed points will exist so long as /ii < 0, and for \i\ > the 
fixed point X^ will exist (and will be unstable) but X|j will not exist. The representative 
bifurcation diagram for such a situation is shown in Fig. [7| When [12 < 0, the period- 1 
fixed point is not observed over the whole parameter range, but the existence of the 
unstable fixed point X^ causes high periodic orbits to occur for fit > 0. On the other 
hand, if ^2 is chosen to be positive, for fi\ < the fixed point ~K.* R will exist but X^ 
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0.4 



x 




0.05 




^1 



Figure 8. A one-parameter bifurcation diagram with /ii as the bifurcation parameter 
and 1x2 = 0.025 . The other parameters tl = 1.5, 5l = 0.3, tr = 0.3 and Sr = 0.3 are 
in R4. 



will not exist and for /i X > both the fixed points will exist. If the parameters are 
chosen such that —(1 + Sr) < tr < 1 + 5r, then is stable and will occur for the 
whole parameter range and is unstable. This is shown in the simulated bifurcation 
diagram of Fig. |H1 

4-3. Conditions for existence of period- 2 orbits 

The conditions for existence of period-2 orbits can be obtained from (j25j) and (|2fij) . First 
we have to find the values of fci, k 2 and k 3 . 
From Section El we have 



A x = A x - 



1 


( bl 


0> 




1) 


[ K 


J 





t l 



p-(-i) 






and 



B = ( ,1 j -(-l)'(--.i I A^J^-Aoh 

l + TR 1 

Sr 

tl - t r + S L - 6 R 
Sr{1 + t l )-S l {1 + t r ) 

d=\l- Mlp*^-!) = (1 + 6 L )(1 + Sr) - T L Tr 



D* = (-l) 2 (adj(I + A 2 ))C 

1 - 1 V 1 

Sr 1 + tr \ 
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Now with the help of above values, we find 
adj(I-Ai)C 1 / 1 



K 



d d \ —5l 

adj(I-Ai)D* 1 



2 



d 

Substituting the parameters, we obtain 




fci 

and from (|2*7)l 



1 



[l + 5 L ){1 + 5 R ) - t l t r 

l + 8 R 
[1 + 8 L )(1 + 8 R ) - t l t r 



(I + 5 L )(l + 5 r ) - t l t r 
Using these expressions, (J25J) and (J2BJ) give 
tr/ii + {l + 5 R )fi 2 



{1 + 8l){1 + 8r)-t l tr 
(1 + 8 L )ni + r L fj, 2 



< (36) 
> (37) 



(1 + 5 l )(1 + 8r) -t l tr 

as the conditions of existence of the period-2 orbit. 

We now show that these are indeed the conditions of the existence of the period-2 
fixed points in the 2D map. 

4-4- Proof of the conditions of existence of period-2 orbits using i2fy) 

Since the system is linear in each side of the border, period two (or a higher period) 
fixed points cannot exist with all points in L or all points in R. However in some regions 
of the parameter space, a period two fixed point may exist with one point in L and one 
point in R. The conditions for existence of period-2 orbit are obtained by solving 

(38) 
(39) 

— •JR u / \ " L 2 I \ u I 

where m\ < and m\* > 0. 

Now solving the above two equations we find 

= W + O^Ml* < o (40) 
1 (1 + 6 L )(1 + 8 R ) - t l t r ■ 
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[1 + 5 L )(Xi + r L/ u 2 



> 



{1 + S L ){1 + 8 R ) - t l t r 

These are the same as those obtained from the theory for n-dimensional systems. 
This implies that if 

T L T R > (l + S L )(l + 5 R ) 

then the period-2 orbit will exist when 

_ T R fM2 (1 + S L ) 

(1 + 5 R ) Hi t l 

and if 



(41) 



(42) 



r L r R < (1 + S L ){1 + 8 R ) 

then the period-2 orbit will exist when 

tr \H ^ (1 + <^l 
Mi 



(43) 



(1 + S R ) H\ t l 

These two conditions are shown in Fig. Efa) as the shaded and unshaded regions 
respectively in the t l —t r space. Now, for particular values of the parameters r L = 0.7, 
5l = 0.3, t r = 0.2 and 5 R = 0.3 the region of existence of period-2 are depicted in the 
ji\—^2 space in Fig. G3^b). 



-(i+8 L )^ 



(i+8 R ) 



(a) 
















(0,0) 







\ (i+5 L ) 



-(1+Sr) 



0.5 



-0.5 
-1 



-0.5 



(b) \ 









0.5 



Figure 9. (a) Regions of existence of period-2 points depending on tl and tr. The 
period-2 points exist in the shaded region if tltr > (1 + 5l)(1 + 6r) and in the 
unshaded region if t^tr < (1 + Sl)(1 + 5r). (b) Regions of existence of period-2 point 
in ^i — fi2 parameter space. The period-2 points exist in the shaded regions for a set 
of particularly chosen tl, tr, Sl and Sr. 



These equations allow one to calculate the parameter range of existence of the 
period-2 orbit, when one of the parameters is varied. To illustrate this, we show the 
bifurcation diagram in Fig. E3J where /ii varied and all the other parameters are kept 
fixed. The stability of the period-2 fixed point is given by the eigenvalues of the 
matrix AiA 2 . For the chosen parameter values, these are Ai = —0.23 + 0.19i and 
A2 = —0.23 — 0.19i This implies that if the period-2 fixed point exists, it will be 
stable. Now from (J3T?)) and (j3*Tj) we find that the period-2 orbit should exist between 
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Hi = 0.01346 and ix\ = 0.1625. Again, since /i 2 < 0, Table [T] predicts that the period-1 
orbit should exist so long as \i\ < 0. The bifurcation diagram clearly shows that the 
above predictions are true. 




-0.1 1 1 1 1 1 1 

-0.05 0.05 0.1 0.15 0.2 



^1 

Figure 10. A representative one-parameter bifurcation diagram with \x\ as the 
bifurcation parameter. The other parameters are: tl = 0.7, 5l = 0.3, = 0.2, 
S R = 0.3, n-2 = -0.025. 



5. Conclusions 

In view of the fact that a large number of systems of practical interest give rise to 
discontinuous maps, it has become necessary to develop a comprehensive theory of 
bifurcations in such systems. In this paper we have taken the first step in that direction 
by extending Feigin's method that was earlier developed in the context of continuous 
maps. We have obtained the conditions of existence of the period-1 and period-2 
fixed points in general n-dimensional discontinuous maps, and have demonstrated the 
application of the derived results in a piecewise linear two dimensional discontinuous 
map. 

It may be noted that the results obtained in this paper speak only of the existence of 
period-1 and period-2 fixed points. In order to obtain information about the occurrence 
of high-periodic orbits and chaos, one has to consider specific dimension of the system, 
which cannot be obtained in the general n- dimensional context. 
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